Critical current of a Josephson junction containing a conical magnet 
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We calculate the critical current of a superconductor/ferromagnetic/superconductor (S/FM/S) Josephson junc- 
tion in which the FM layer has a conical magnetic structure composed of an in-plane rotating antiferromagnetic 
phase and an out-of-plane ferromagnetic component. In view of the realistic electronic properties and magnetic 
structures that can be formed when conical magnets such as Ho are grown with a polycrystalline structure in 
thin-film form by methods such as direct current sputtering and evaporation, we have modeled this situation in 
the dirty limit with a large magnetic coherence length (£/). This means that the electron mean free path is much 
smaller than the normalized spiral length X/2n which in turn is much smaller than (with A as the length a 
complete spiral makes along the growth direction of the FM). In this physically reasonable limit we have em- 
ployed the linearized Usadel equations: we find that the triplet correlations are short ranged and manifested in 
the critical current as a rapid oscillation on the scale of \/2tt. These rapid oscillations in the critical current are 
superimposed on a slower oscillation which is related to the singlet correlations. Both oscillations decay on the 
scale of We derive an analytical solution and also describe a computational method for obtaining the critical 
current as a function of the conical magnetic layer thickness. 



I. INTRODUCTION 

The interaction of singlet-type superconductors (S) with 
ferromagnetic materials in S/FM hybrid systems is a field of 
extensive and ongoing research (see Refs. 1,2,3 and references 
therein). In proximity, the interaction of these competing elec- 
tron orders is characterized by an oscillating component in 
the Cooper pair wave function which leads to a number of 
interesting phenomena: the critical superconducting tempera- 
ture T c dependence of S/FM bilayers on FM layer thickness 
rfj j 4 i 5 i 6 i 7 dependence of T c on the orientation of FM layers 
in FM'/S/FM" spin valvea 8 ' 9 ' 10 ' 11 ' 12 ' 13 and S/FM'/FM" mul- 
tilayers, and finally the realization of 7r coupling in S/FM/S 
Josephson junctions . 14 ' 15 ' 16 ' 17 ' 18 

The standard analysis of the S/FM systems has mostly as- 
sumed that the FM is homogeneous and collinear, in which 
case only the singlet superconducting correlation appears in 
the theory. Extending this standard approach, theory strongly 
indicates that if the FM is inhomogeneous and noncollinear, 
the longer-ranged triplet superconducting correlations should 
then emerge at the S/FM interface. 1,2,3 These triplet correla- 
tions should then be insensitive to the exchange field of the 
FM material and as such their proximity range is expected to 
be similar to that of singlet pairs in a superconductor/normal 
metal system. 

Inhomogeneous magnetization exists in a range of mate- 
rial systems, which can be classified into three categories: (1) 
magnetic domain walls; (2) ferromagnetic multilayers such 
as when FM layers are decoupled via a nonmagnetic (NM) 
spacer to form spin-active devices; and (3) the intrinsically 
inhomogeneous and noncollinear magnetic materials. 

Domain walls were one of the first magnetically inho- 
mogeneous systems to be combined with superconductivity. 
Although experimental studies of such systems are notori- 



ously challenging because of the need to control the mag- 
netism at the nanometer scale, results and analysis have in- 
dicated that domain walls are favorable nucleation sites for 
superconductivity ! 19 ' 20,21,22 ' 23 Theoretically, the emergence of 
triplet components in junctions containing a single domain 
wall and or a multidomain ferromagnet (MDFM) have been 
extensively analyzed. 24,25 Recently, large area S/MDFM/S 
junctions have been fabricated. 26 In this type of junction, the 
amplitude of the critical current is expected to decay expo- 
nentially with FM layer thickness. If singlet-type electron 
pairs are scattered into triplet ones at the domain-wall regions, 
it is expected that for a critical thickness of the MDFM the 
triplet correlations will dominate over the singlet ones leading 
to slower decay in the critical current with the MDFM thick- 
ness. So far, evidence of a crossover from singlet to triplet- 
dominated transport in these types of systems is nonexistent. 

The second category, the ferromagnetic multilayers, have 
been combined with superconductors in S/FM'/FM'/S junc- 
tion form junction form although most studies have been the- 
oretical up to now 27,28,29,30 with only a few experiments show- 
ing how the Josephson ground state is affected by the ori- 
entation of the FM layers. 31,32 The majority of experimen- 
tal studies have focused on how a superconducting layer is 
modified by the relative orientation of the FMs. In these sys- 
tems, however, the triplet superconducting components that 
exist when the FMs are noncollinear only transmit informa- 
tion about the direction of the magnetic layers. 33,34 To observe 
a longer-ranged spin triplet proximity effect, it is thought that 
the Josephson junction must contain three more FMs ; 30 ' 35 with 
each offset from the other by an angle 6 ^ [0, ir] (with tt as the 
antiparallel configuration). In principle, the angle 8 and thus 
the triplet components could be controlled by the application 
of an external magnetic field. Unfortunately, the implementa- 
tion of a large enough change in the angle with an applied 
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magnetic field is very difficult to realize without strongly sup- 
pressing the superconductivity. 

The third category, the intrinsically noncollinear magnets, 
is potentially one of the simplest systems to combine with a 
superconductor to experimentally study triplet correlations. 36 
Recently, 37 interferometer measurements of superconducting 
Al coupled to the rare-earth metal Ho have been made. In 
these Al/Ho/Al junctions, superconducting phase periodic 
conductance oscillations were observed indicating the pres- 
ence of a longer-ranged proximity effect when interpreted in 
the limit of a small coherence length in the Ho relative to the 
length of a complete spiral A. 38 It is understood that the triplet 
correlations were generated at the Al/Ho interface due to a 
rotating magnetization present there and sustained by a con- 
tinuous magnetic spiral throughout the length of the Ho. A 
similar explanation 39 ' 40 ' 41 ' 42 was given for a long-ranged prox- 
imity effect observed in the half-metal CrC»2. 43 In this system, 
the triplet current was shown to be insensitive to the strong 
polarization of the half metal. Spin mixing at the interface is 
currently the best explanation for the triplet proximity effect 
observed although a better understanding of the interfaces that 
can exist in these types of material systems is needed to ver- 
ify this explanation. For Ho, it is well known that growing 
it in thin-film form with a magnetic spiral at the interface is 
difficult to achieve. This again highlights a need to improve 
our understanding of the likely properties and structures that 
can arise at the interface of noncollinear magnets, such as Ho, 
with superconducting materials. 

The magnetic structur e) 44 ' 45 ' 46 and electronic/thermal 
properties 45 of the rare-earth Ho are well known. Its magnetic 
structure has been characterized in bulk, single crystal, and 
thin-film forms by neutron diffraction, x-ray diffraction, 
and vibrating sample magnetometery. In thin-film form, the 
quality of the conical magnetic structure is poorly understood 
although it is well known that the growth method and growth 
conditions, crystal forms, and interfacing materials affect the 
ordering range of the magnetic structure^ 4 ^ 5 - 

Long-ranged magnetic ordering in Ho requires a coherent 
crystal structure in which the c axis is the screw axis with 
the moments in the basal plane configured into a distorted 
helix parallel to the c axis. The quality of the Ho (e.g., im- 
purity content and roughness) and the strain at the NM/Ho 
interface are both important factors in determining the scale 
of magnetic ordering; for example, substantial intermixing 
at the NM/FM interface may disturb the growth in the he- 
lix which may affect, smear out, or even destroy any triplet 
correlations. Neutron-diffraction studies on epitaxial (inter- 
facially strained) Nb/Ho bilayer films grown by dc magnetron 
sputtering 47 at high temperature suggest the presence of an in- 
plane spiral (antiferromagnetic part) but no out-of-plane pitch 
(ferromagnetic part) was detected even down to very low tem- 
peratures T ~ 1 K. This implies that the strain at the Nb/Ho 
interface is suppressing the ferromagnetic component. Fur- 
ther studies on polycrystalline Nb/Ho/Nb trilayer films have 
also been made. 48 In these films strain at the Nb/Ho interface 
is lower and from Josephson-junction-type measurements a 
weakly conical magnetic structure was confirmed from field- 
dependent measurements of the junction's critical current as a 
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FIG. 1: (Color online) (a) An illustration of two singlet-type super- 
conductors (S) sandwiching a ferromagnet with a conical magnetic 
structure (Ho). The two magnetic phases of Ho, (b) the in-plane an- 
tiferromagnetic spiral phase of Ho below its Neel temperature O n , 
and (c) the conical magnet phase of Ho below its Curie temperature 
<d c : the antiferromagnetic component I r rotating in the {x, y} plane; 
the ferromagnetic component I z pitched towards the z axis; and the 
resultant magnetization vector I rotating on the surface of a cone. In 
the limit considered in this paper, the ferromagnetic coherence length 
is much larger than X/2n. 



function of the Ho spacer layer thickness. 

Altogether, the above review shows that although Ho can 
be grown on top of thick Nb leads with a conical magnetic 
structure, strain at the interface does weaken the ferromag- 
netic component possibly implying a weaker magnetism. The 
motivation behind this work is to complement the currently 
available theory on S/FM/S junctions with a conical FM weak 
link by considering the physically reasonable situation (see 
Sec.|nl| in which the conical ferromagnetic coherence length 
£f is much longer than the normalized spiral length X/2tt. 
Because both and \/2ir are much larger than the electron 
mean free path I (dirty limit), the S/FM/S junction can be de- 
scribed within the framework of the linearized Usadel equa- 
tions. 

This paper is organized as follows: Sec. |II]reviews the mag- 
netic and electronic properties of thin-film Ho and outlines 
the important physical properties of the situation being ana- 
lyzed in this paper; in Sees. [IIT1and [TVl the general theoretical 
framework in which we model a Josephson junction contain- 
ing a conical magnet is described with analytical solutions to 
obtain the Josephson current explained; in Sec. [VJ we present 
a computational method for calculating the Josephson current, 
which is particularly useful to experimentalists. 



II. MAGNETIC STRUCTURE AND ELECTRONIC 
PROPERTIES OF THIN-FILM HOLMIUM 

Consider the general structure of a conical magnet which 
consists of a rotating in-plane magnetization and a con- 
stant out-of-plane magnetization [see illustrations in Figs. 
02 a )-{n c )]- The in-plane component is effectively an antifer- 
romagnetic (AFM) state which orders itself at the Neel tem- 
perature 8„, while the out-of-plane component can be con- 
sidered to be a ferromagnetic phase which orders itself at the 
Curie temperature of the material. Thus, the strength or ex- 
change interaction energy / of the ferromagnetic part is re- 
lated to the Curie temperature: / ~ fcg6 c . The in-plane com- 



3 



ponent completes a full rotation in a distance of A along the z 
axis, which implies the distance along z on which the in-plane 
component rotates in 1 rad is X/2tt. From now on, X/2tt will 
be referred to as the normalized spiral length. 

In the analysis that follows this section we shall assume that 
the electron mean-free path i is smaller than both the coher- 
ence length of Ho and the normalized spiral length X/2n. 
For Ho in thin-film form, this limiting situation is justified for 
the case when it is sputter deposited and polycrystalline. 

Polycrystalline thin films of Ho have a large residual resis- 
tivity po in the (6— 12) x 10~ 7 Om range (see Refs. — and^and 
references therein). A rough estimate of the electron mean 
free path for the conduction electrons around 4 K using the 
relation £ — vftti/ p^ne 2 , where vp, the Fermi velocity is 
1.6xl0 6 m/sr^ m is the mass of the electron, and n is the 
number density of free electrons, gives a (0.5—1 .0) nm range, 
which is smaller than both £f ~ (6—7) nm and X/2tt ~ 1.1 
nm4£ Even in single-crystal form, the resistivity of thin-film 
Ho is large and around 6x 10 -7 Sim with an electron mean- 
free path of ~ 1.0 nm in the c-axis orientation. 



III. GENERAL SOLUTION IN A CONICAL 
FERROMAGNET 

Let us consider a conical ferromagnet FM, the axis of which 
coincides with the z axis. The magnetization vector and hence 
the exchange field I = (I x ,I y ,I z ) has a constant axial compo- 
nent I z and a radial component I r , which rotates in the {x, y} 
plane with wave vector Q = 2ir/X as we move along the z 
direction (see Fig. [TJ. The x and y components of I are there- 
fore given by 

Ix{z) — I r cos{Qz) and I y (z) — I r sm(Qz). (1) 

In this section we find the solutions for the anomalous Green's 
function in the conical FM, considering in particular the case 
of large Q. We assume that the dirty limit is fulfilled, which 
means that the FM coherence length (£/ = \JhDj /|I| with 
D f being the diffusivity of the FM) and the normalized spi- 
ral length X/2tt are both much larger than the electron mean 
free path, i.e. A/27T ^S> £■ If it is also assumed that the 
anomalous function is sufficiently small in the FM (which is 
the case if the S/FM interfacial resistance is large enough), 1 
we can use the linearized Usadel equation^ 



TiD 



d 2 F{z) 
dz 2 



-2n\u\F{z)-sea{u)i\l{zy,F{z)] + = 0, (2) 



substitute expression © into Eq. (O we obtain a set of equa- 
tions for the four components of the anomalous function F, 

~hD f ^^-huf (z)-i[I(z)-f(z)]=Q, (4) 



\ hD f d2fx d V z 2 {Z) -*of*,vA*)-iI*,vA*)M*) = 0- (5) 

Since the symmetric properties of the Usadel equations with 
respect to w are trivial, only the case of u> > is treated from 
now; we already omitted sgn(w) and used w instead of \u>\ 
in Eqs. (01 and ((5]). After putting expressions (Q]) into Eqs. 
(0]l and Q, they can be simplified with the substitution f± = 
f x ±if y , which yields 



dz 2 



-2klf -i(2kif z + k' r (f + e^ z + f^ z ))=0, (6) 



dz 2 



2k 2 J z - 2ik 2 f Q = 0. 



(7) 



(8) 



The quantities fc w = u> J D / and k z ^ r — I z r /HD f were 
also introduced at this step. 

By searching a solution for Eqs. ©—([8]) in the form of 



k,z=A Q<z e Kz and /± = A ± e" z e 



Kz„±iQz 



(9) 



with A± = A' + ± A'_, we obtain a set of algebraic equations 
for the amplitudes Aq, A z , A' + and A'_: 

{K 2 - 2kl)A - 2ik 2 z A z - 2ik 2 A' + = 0, (10) 



2ik 2 A + (K 2 - 2kl)A z = 0, 



(11) 



2ikt.Au + (K 2 - Q 2 - 2kl)A' + + 2iKQA'_ = 0, (12) 



2iKQA' + + (K 2 - Q 2 - 2kl)A'_ = 0. 



(13) 



A nontrivial solution only exists for the amplitudes if the de- 
terminant of the system is zero; this gives a fourth-order equa- 
tion for K 2 , 



[{K 2 - 2kl) 2 + 4kt] [(K 2 - Q 2 - 2kl) 2 + 4K 2 Q 21 



where the anomalous Green's function 



F = / l + f 



(3) 



is a matrix in spin space with f = (f x , f y , f z ), and a — 
(a x ,a y ,a z ) is a vector containing the Pauli matrices. The 
component /o is an even, while f x , f y and f z are odd func- 
tions of the frequency u>. The Matsubara frequencies are given 
by uj = (2n + l)nk B T/h with n = 0,±1,±2, ... at temper- 
ature T and [a, b] + = ab + ba is the anticommutator. If we 



+4k 4 r (K 2 - 2k 2 J(K 2 -Q 2 - 2kl) = 0, 



(14) 



which is equivalent to the similar equation obtained by Volkov 
et alM- In their paper, they considered the limit in which 
k r ,k z S> k u , Q, whereas we take the limit of Q ^> k r ,k z ,ku;. 
This limit seems to be appropriate in the case of Ho, which 
has a conical magnetic structure with A ~ 6 nm and there- 
fore Q ~ I nm -1 . The exchange energies I r and I z can be 
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estimated from the AFM and FM ordering temperatures; as- 
suming a typical diffusivity Dj sw 5x 10~ 4 m 2 s _1 and a tem- 
perature T w 4 K, we obtain fe r ~ 0.2 nm" 1 and fc 2 , fc w ~ 
0.05 nm" 1 , which are all much smaller than Q. 

In order to make the approximations more transparent, we 
introduce the dimensionless quantities 

K 2 k r z u 

A = and = ' ' (15) 

Q 2 <y 

with e r z w <C 1. Without assuming anything about the rela- 
tive values of these small numbers, we take the case in which 
their respective leading terms are on the same order of mag- 
nitude; our results are therefore applicable to the general case 
and the particular cases can be obtained by taking appropriate 
limits. It turns out that the leading terms in the small quanti- 
ties e, i2 W are on the orders of e 2 , e 2 and e u , hence we assume 
for the approximations that e 2 ~ e 2 ~ e w . Two roots of A 
are on the order of 1 ; if we neglect every term smaller than e 2 
from Eq. ( fT4b . we obtain 

(1 + A) 2 = ^(1 + A + 2A 2 ) + ^(1-A). (16) 
A A 

Since the terms on the right side are <C 1, the left side has 
to be small, which is only possible if A w —1. In this case 
we can substitute A = — 1 on the right side, hence Eq. ( fTBT l 
reduces to 

(l + A) 2 + 8e w + 8e 2 = 0, (17) 

which gives A = — 1 ± i^J?>e u + 8e 2 . If we only keep the 
roots of K for which Kc(K) > and still neglect the terms 
smaller than e 2 , we obtain 

K h2 = ±iQ + Q^2e^ + 2e 2 (18) 

for the first two eigenvalues K. These correspond to rapidly 
oscillating solutions (together with the rotation of the magne- 
tization vector), which decay much more slowly in the neg- 
ative direction. Since K only appears as K 2 in Eq. ( TT4"1 >. 
the roots with Re(A") < can all be paired up with their 
respective opposites and give the same solutions decaying in 
the positive direction. Expression (fT8l without the term 2e 2 is 
equivalent to the result obtained by Bergeret et al. 2 for a spiral 
ferromagnet (I z = 0, hence e 2 = 0). 

The two remaining roots for A are on the order of e 2 ; in this 
case we can treat A as being small and hence neglect larger 
powers of it. However, we must keep terms up to the order of 
ef, in Eq. ( TBI ) to obtain the quadratic equation 

A 2 - 4(e 2 + e w )A + 4(e 2 + e 2 + 2e w e 2 ) = 0, (19) 

which yields A = 2e UJ + 2e 2 ± 2^/e 4 . - e 2 . Note that the 
leading terms are indeed on the orders of e 2 , e z and e w , as 
stated above. Two more eigenvalues K with Rc(K) > are 
obtained, 

K 3A = Q^2e^ + 2e 2 ±2^4-e 2 . (20) 



The behaviour of these solutions depends on the relative val- 
ues of e z and e 2 and now we can consider the two particular 
cases. If e 2 > e 2 , the roots K 3A are complex conjugates and 
the solutions (l20l describe a slowly decaying oscillation in the 
negative direction. If e z < e 2 , the roots K3 4 are real, which 
means that F decays exponentially without oscillations. This 
case corresponds to almost in-plane magnetization and con- 
tains the limit of th e spiral fer romagnet; expression (f20b re- 
duces to K 3 = Qy / 2e UJ + 4e 2 and K± = Q^/2eZ if e 2 = 0. 
The solution corresponding to K4 has zero amplitude in any 
S/FM system, while K 3 coincides with the value obtained by 
Bergeret et al— 

After determining the eigenvalues K we calculate the cor- 
responding eigenvectors, i.e. the relative amplitudes of the 
different components /o, f z and f± in each solution. Since the 
roots -Kj.,2 given by Eq. ( fT8l appear as a direct consequence of 
the rotation of the magnetization vector, we expect the com- 
ponents f± to dominate in the corresponding solutions, and 
hence we choose A' 1+ = A' 2+ = 1 [A' 1+ and A' 2+ are the A' + 
amplitudes appearing in Eqs. (fTOb — (fT3b for the solutions cor- 
responding to the roots Ki and K 2 , respectively]. Equation 
( fTTT ) shows that A z <C Aq in these cases, while Aq <C A' + = 1 
according to Eq. ( fTOb . It is valid therefore to take A z w 0, 
then use Eqs. ( fTOb and (TT~3T > together with Eq. ( TT~8T > to obtain 
A[_ = —1, A 2 _ = 1 and Aw = A 2 q = —2ie r in the leading 
approximation. 

The solutions corresponding to the other two roots K3 4 
predominantly consist of the components /□ and f z , there- 
fore we choose A 30 — A 40 = 1. Equations dT2b and ( fT~3b 
show that A'_ <C A' + <C Aq in these cases, which implies that 
A'_ ~ 0. Keeping this in mind, we can apply Eq. ( fT2b to get 
A' 3+ = A' i+ = -2ie r , and Eq. (Qj) with Eq. (|20j to obtain 

A 3z = f== and A iz = f==- (21) 

We can again consider the two cases: if e z > e 2 , A 3z and A4 Z 
are complex numbers with unit modulus, and A± z — —A* 3z . 
In particular, Eq. (fJTJ reduces to 

2 2 
A Zz K,\ + i— and A Az w -1 + i— (22) 

in the limit of e 2 3> e 2 . If e 2 < e 2 , ^3 2 and ^4 2 are purely 
imaginary numbers with |^43 2 ||^44 2 | = 1; they are being ap- 
proximated as 

A 3z * 1 ^ and A 42 «^ (23) 
2a 2 e 2 

if e 2 -C e 2 . In the limiting case of the spiral ferromagnet 
(e 2 — > 0), A 32 — * and |^4 2 | — + 00. The latter means that 
the solution for F corresponding to the root K4 has only f z 
component (and no singlet /o component); its amplitude is 
therefore zero, as already mentioned above. 

If we take the eigenvalues K with Ke(K) < 0, the corre- 
sponding eigenvectors are similar to those corresponding to 
their opposites; Eqs. (fTOb — <fl~3T> show that the amplitudes A , 
A z and A' + remain the same if we multiply K by (—1), while 
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A'_ changes sign. According to A± = A' + ± A'_, this means 
that A + and A_ are exchanged. Keeping this in mind, we 
can write down the general solution for the linearized Usadel 



equation (O in a conical ferromagnet. If the FM occupies a 
region of thickness d in the z direction (more specifically, the 
range < z < d), the general solution can be written as 



F = ^\B n e- K ^{A n0 l+A nz a z + -Y,A n Te ± ^\a x Tia y ^ 

I l ± z ± 



where the eigenvalues K n and the relative amplitudes A n $, 
A nz and A n ± — A' n+ ± A' n _ are given by the expressions ob- 
tained in this section. Note that the amplitudes A n+ and A n - 
are exchanged in the terms corresponding to the solutions with 
Ke(K ) < 0, as mentioned above. The amplitudes B n and 
C n are determined by the boundary conditions at z = and 
z = d; these are discussed in Sec. HVl 
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IV. JOSEPHSON CURRENT IN THE S/FM/S JUNCTION 

If the regions with z < and z > d on the two sides of the 
FM are occupied by two identical half-infinite superconduc- 
tors S, we obtain a S/FM/S junction. The axis z of the conical 
FM is perpendicular to the S/FM interfaces (see Fig. [TJ. The 
two superconductors have a phase difference of $ with respect 
to each other, so the bulk pairing potentials in the left and the 
right S are given by Ae~ 4<1>/2 and Ae 1 */ 2 (A € M). The nor- 
mal and the anomalous Green functions are Gl,r — G s l and 
F l ,r = FsleT^/ 2 in the bulk of the left and right super- 
conductors, respectively, where G s = huj/ ^Jh 2 uj 2 + A 2 and 
F s = A/ \Jh 2 L0 2 + A 2 . The normal state conductivities of the 
S and the FM are a s and <j/, while the interfacial resistance 
per unit area between the S and the FM is denoted by R. We 
introduce the dimensionless quantities 



FIG. 2: (Color online) Typical I c R n dependence on the thickness d 
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f 



(25) 



where £ s = y/hD s /2TTksT is the superconducting (quasipar- 
ticle) coherence length of the S with D s being its diffusivity. If 
the interfacial resistance is large enough, i.e. 7^ ^ max(l, 7), 
we can use rigid boundary conditions at the S/FM interface; 1 
we assume that the pairing potential and hence the Green's 
functions are the same at the interface as in the bulk material. 
Furthermore, because of 76 ^> 1 the anomalous function F 
is sufficiently small in the FM, which verifies using the lin- 
earized Usadel equations in the FM (see Sec. HlB . 

Assuming that 7^ ^> max(l, 7) is true, the rigid boundary 
conditions are 1 



F L = G s F(0) - 7 B C/ 
at the left side of the FM (z = 0) and 

F R = G s F{d)+ lB £ f 



dF(0) 
dz 



dF{d) 
dz 



(26) 



(27) 



of the FM layer (A 

xl0~ 4 mV\ D s 

n6 



= 2.2x10" 
2.5 xlO 



J, Q 



-- 9 xl(f m -i , 
a f = 4 x 10 6 



Dj = 6 
(Jim)" 1 , 

cj s = 6 xlO 6 (finiT 1 ): (left) in the limit of e z > t 2 r (I r /k B = 
h/k B = 100 K) and (right) in the limit of e z < e 2 r (I r /k B = 130 
K, I z /kB < 4 K) for two different values of the interfacial resistance 
R. 



at the right side of the FM (z = d). If the FM layer is thick 
enough (d ^s> £/), the terms containing eT K ^ d ~ z ) can be 
neglected from the general solution (l24l > near 2 = 0, while the 
terms containing e~ K,lZ can be neglected near z = d. In this 
case we can take the components f , f z and f± of Eq. ( |24| ) at 
z = d and obtain equations for the amplitudes C n : 



— e 1 * 72 = 



n=l 



C n A nQ (l+TK n ), 



Q = Y d C n A nz {l + TK n ), 



= Y,C n A n± (l + r{K n ±iQ)), 

n=l 



(28) 



(29) 



(30) 



where the notation T = 7b^/ /G s is used. Similar equations 
hold for the amplitudes B n at z = 0, the only difference is 
that the sign of the phase <£>/2 in the first term of Eq. (|28T > is 
negative. Substituting K n and A n ± into Eq. ( f30b and using 
I K3.4 1 < Q yields 

(d, 2 - ie r (C 3 + C 4 )) + T(±Qe r (C 3 + C 4 ) + k C xa ) = 0, 

(31) 

where k = Qy/2e LU + 2e 2 . Even though TQ > 1, TQe r < 
1 for realistic values of the interfacial resistance; it follows 
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from Eq. d3TT l that C\,C2 <C C 3 ,C 4 , therefore the terms 
containing C\ and C 2 can be neglected from Eqs. d28l > and 
Those equations with A 3 q = A 4 q = 1 take the form of 



A_ 



!$/2 = (C 3 + C 4 ) + T(K 3 C 3 + K±C 4 ), 



(32) 



- (A 3z C 3 + A 4z C 4 ) + T(A 3z K 3 C 3 + A 4z K 4 C 4 ). (33) 

In the following we work in the limit of e z ^> e 2 , so we take 
A 3z and A 4z given by Eq. d22l) to get 



Q 



3,4 



Ae'*/ 2 / .e 2 .. 

2Mi + ^3, 4 rH T \ z 1 



(34) 



where the expression (|20b for the roots ^3 4 can be approx- 
imated as if 3. 4 « Q^/2e u) + 2e 2 ± 2ie z in this case. Sub- 
stituting C 3 and C 4 into Eq. (|3TT l gives the remaining two 
amplitudes, 



Ae' $ / 2 £r (»Trg) / l-ie 2 /e z 1 + ie 2 /e, 
2fiw(l + fc r) ^ 1 + K 3 T 1 + K 4 T 



(35) 

The term i can be neglected from the numerator, because 
TQ ^> 1. Furthermore, since K 3 and K 4 are complex conju- 
gates if e z > e 2 , the sum in Eq. d35l l can be simplified to 



^ Ae^/ 2 e r TQ Rc (l-ie 2 r /e z 



Huj(1 + k T) 



i + K 3 r 



(36) 



The results for B 3 ,4 and B12 are the same as those given by 
Eqs. d34l > and (f36l > with the only difference being a minus sign 
before z<E>/2 in the phase factor of A. 

The Josephson current in the S/FM/S junction of area A is 
given by 35 



TTUfA . 



■ B T^lm\^(F*{z)a, 



UJ>0 



dF(z) 
dz 



(37) 



which can be evaluated in the range < z < d. We can 
substitute the solution d24l i into Eq. ( f37T > and take z = d; 
most terms do not give any imaginary contribution to the trace, 
hence we obtain 



/= Z^ fcsr Vlm(S 1+ S 2 ), 

P ' 



(38) 



S 2 = -J2(j2 B n A nT^ Kd J2 CnA ^ Kn±l ^ 



4 



4 



C* A 

/ , ^n-^n^f / j 



B n A n ±(K n ±iQ)e 



-K n d 



(40) 



The first term S\ mainly contains the solutions corresponding 
to the eigenvalues K 3 ,4, while the second term S2 is mainly 
contributed by the solutions corresponding to Ki 2 . By using 
the values of A n o, A nz and A n ± in the limit e z ^S> e 2 , the 
expressions for Si and S% become 



.V, = 4K 3 e- K3d [1+if] (B* 4 C 3 - C* 4 B 3 



+ 4K 4 e 



-K d d 



l-i^)(B* 3 C 4 -C* 3 B 4 ), (41) 

€2 



S 2 = 4fc [e- Xld (C 2 *B 1 -B* 2 C X ) 



,-K 2 d 



(C{B 2 



■B* X C 2 )] 
(42) 

in the main approximation. Here we neglected the terms con- 
taining the small amplitudes A 4 o, A 2 q, A 3 ±, A 4 ± and used the 
fact that \ A 3z \ = \ A 4z \ = 1 if e z > e 2 . By taking K 2 = K{* 
and K 4 = K 3 into account, then putting the above obtained 
expressions for B„ and C n into Eqs. (l4"TT i and d42b we obtain 



A 2 

Si = 4isin($) — —Re 



K 3 e 



-K,d 



(1 + K 3 vy 



(43) 



16i sin($)- 



Re 



1 - ie 2 r /e z 
l + K 3 T 



e 2 Q 2 r 2 fc e- fc ° rf 

(i + k vy 



(44) 



0,z n—1 n—1 
4 4 

- T, C n A n(o^)T, B -Mo, z )K n e- K - d ), (39) 



n=l 



n=l 



The Josephson current through the S/FM/S junction therefore 
obeys the formula I = I c sin(<l>), where the critical current I c 
is given by 
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I c Rn = 4nk B T 



d + 2£ /7B 



w>0 



h 2 u 2 



Re 



K 3 e- 



-K 3 



(i + ^ 3 r) 



r 



Re 



1 - itr/ e z 

i + K a r 



(i + hvf 



- ka d 



■ cos(Qd) 



J45} 



with i?„ = (d + 2^/7b)/ct/A being the normal state resis- 
tance of the junction. 

The dependence of I c R n on the thickness d predicted by 
Eq. d45T > is plotted in Fig. |2](left) for two values of the interfa- 
cial resistance R. Both curves show a small, rapid oscillation 
superimposed on a large, slow oscillation; they both decay on 
the scale of the slow oscillation. Comparison between the two 
curves demonstrates that an increase in R reduces the current, 
but makes the rapid oscillations relatively more pronounced. 

The first term in Eq. d43T > gives the slow oscillation, which 
is mainly due to the "short-range" singlet and triplet compo- 
nents of the anomalous Green's function F (i.e. the singlet 
component fo and the triplet component f z with zero projec- 
tion on the z axis). Conversely, the second term in Eq. d45l l 
corresponds to the rapid oscillation, which is related to the 
"long-range" triplet components (i.e. the triplet components 
f± with projection ±1 on the z axis). Note that in our case 
the terms "short-range" and "long-range" do not mean any 
difference in the respective decaying lengths; they are only 
defined like this to be consistent with the notions used in other 
papers , i 

Unlike the slow oscillation which is also present in a sys- 
tem with a homogeneous FM, the rapid oscillation appears 
as a direct consequence of the inhomogeneous magnetization. 
This is shown clearly by the coincidence of its oscillation pe- 
riod and the magnetic spiral wavelength A. The magnetization 
changes quickly with respect to the FM coherence length £/, 
which explains why the amplitude of the rapid oscillation is 
small compared to that of the slow oscillation. 



By taking the limit of I r — ► (and hence e r — ► 0), we 
recover an S/FM/S junction with a homogeneous FM of ex- 
change energy I z . In this limit, Eq. (|45T > reduces to 



I c Rn = 4irk B T 



d + 2£/7B 



A 2 

E^ Re 

u>0 



K 3 e- K3d 



(i + K 3 vf 



with the root K 3 taking the form of 



K 3 = y/2k*+2ikl 



(46) 



(47) 



This is the standard formula for I c R n in a Josephson junction 
with a homogeneous FM weak link^ 

The same result is obtained in the limit of Q — > oo because 
e r = k 2 /Q 2 — ► in this case. The amplitude of the rapid 
oscillation vanishes as Q~ 2 , and hence we recover Eq. d46l ). 
This is physically understandable; as the FM coherence length 
becomes very much larger than the characteristic spiral wave- 
length, the radial magnetization "averages out" on the scale of 
£f, which means that the situation is equivalent to I r — ► 0. 

Now we can return to Eqs. d32l and d33l > and take the op- 
posite limit, i.e. where e z <C e 2 .. In this case we use A 3z and 
At± z given by Eq. ( f23l ) to obtain different values for the am- 
plitudes B n and C n . However, the expression ( l38l l with the 
same terms Si and S% still holds for the Josephson current. 
After substituting the new values of B n and C n into Eq. d38b 
and taking approximations valid in the given limit, we recover 
I = I c sin(<l>) and obtain 



J 



I c R n = 4:irk B T 



A 2 

Y — 

UJ>0 



K 3 e 



-K 3 d 



(i + K 3 ry 



K 4 e 



-Kid 



K (i 
i — 



KiT) 2 



4e 2 Q 2 r 2 fc e- 



■ ko d 



(l + k Q T) 2 (l + K 3 T) 2 



cos(Qd) 



(48) 



for the critical current. The roots X3.4 given by Eq. (l20~b can 
be approximated as K 3 



Q y /2e UJ + 4e 2 r and K 4 « 
if e z < e 2 . 

The I c R n dependence on d as given by Eq. d48b is repre- 
sented in Fig. 12 (right). The rapid oscillation is similar as in 
the limit of e 2 ^> e 2 , but the slow oscillation is absent; the 
other component of I c decays exponentially without oscilla- 
tion. The rapid oscillation is still related to the "long-range" 
triplet components, whereas the exponential decay is to the 
"short-range" singlet and triplet components. 

In the case between the two limits (where e z ~ e 2 ), both ex- 
pressions (|43T > and ( |48T > are applicable, but they are not as ac- 



curate as when they are used in their respective limiting cases. 
Since the approximations leading to Eq. (|43T > are less sensitive 
than those required for Eq. d48l i. the former is preferred to be 
used in such a case. 



V. COMPUTATIONAL METHOD FOR CALCULATING 
THE JOSEPHSON CURRENT 

In this section we describe an alternative method for ob- 
taining the Josephson current in the S/FM/S junction; it re- 
quires computational power and does not yield an analytical 
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formula, but is exact within the framework of the linearized 
Usadel equations. The basic steps are the same as in the pre- 
vious sections: we first solve the linearized Usadel equation 
(O with the boundary conditions [Eqs. (|26| > and d27lll, then 
evaluate Eq. d37l i at a suitable location. 
Let us introduce the formal vector 



F(*) = (/<,,/« 



fyifz: 



dfo dfx 
dz ' dz 



df y df z 
dz ' dz 



(49) 



containing the components of F and their respective deriva- 
tives with respect to z. We denote the value of this vector 
at the right side of the left S (at z = 0) and at the left side 
of the right S (at z = d). These consist of the components 
of Fl and Fr, respectively. By using the method described in 
the Appendix, we can relate F^ and F# through a matrix type 
equation. Since we know the first four components of both F^ 
and Ffl, we can use this equation to obtain the remaining four 
components (the derivatives). 

The Josephson current is evaluated with Eq. (|37T i in the S 
side of the left S/FM interface; since we calculate the current 
in the S, we must substitute a s instead of 07 in Eq. J37t . By 
using f X L = f y L = fzL = the formula simplifies to 



I = 



2na.A 



fc B T£lm(/ 0i ^Y (50) 



where foL — Ae~**/ 2 /-\/A 2 + h 2 uj 2 and dfoL/dz is calcu- 
lated by the method described in the Appendix. By setting the 
phase difference to $ = ir/2 we obtain 



I c R n = 2nk B T 



a s {d + 2i n B ) 
ea-f 



w>0 



iOL- 



dfoL 

dz 



(51) 



Note that the values of / z, and dfoL /dz depend on the phase 
difference <£>, as well as on other parameters describing the 
junction. 

Evaluating rf/oi /dz requires inverting matrices and taking 
matrix exponentials, therefore this method does not give an 
analytical formula like expressions ( |45T > and d48l i. On the other 
hand, it gives the right result in the more general case, even if 
Q is not large [in which case neither Eq. d45b . nor Eq. (08]) is 
applicable]. This method can also be used to check analytical 
results; in our case it seems that within their respective ranges, 
Eqs. (|45]> and d48l ) show good coincidence with the results 
obtained by the computational method: see Fig. [3] 



VI. SUMMARY 

We have calculated the Josephson current in a superconduc- 
tor/ferromagnetic/superconductorjunction in which the ferro- 
magnet has a conical magnetic structure. In view of the realis- 
tic interfaces that can exist between thin-film superconductors 
(Al, Pb, and Nb being the elements typically used) and thin 
films of conical magnets such as Ho, we have extended the 
problem to a regime in which the ferromagnetic coherence 
length is long compared to the electron mean-free path and 
the normalized spiral length of the magnetic spiral. 





FIG. 3: (Color online) Comparison of the analytical (solid curve) and 
computational (dashed curve) results for I c R n in the function of the 
thickness d: (left) in the limit of e z 3> and (right) in the limit of 
e z <C tr- The parameters used are identical to those in Fig. |2]with 
R = 3.0xl(T 15 ftm 2 . 



From the materials point of view, the dirty-limit model we 
present is physically reasonable and most applicable for when 
the Ho thin film is polycrystalline. The electron mean-free 
path of thin film Ho is < 1.0 nm and its normalized spiral 
length is ~ 1.1 nm, whereas the coherence length in such films 
has recently been determined to be in the 6—7 nm range. 48 

In this new situation, we have shown that the Josephson cur- 
rent is highly sensitive to the length of the conical ferromagnet 
A with the current containing a rapidly oscillating component 
in the function of the total conical magnetic thickness. These 
rapid oscillations are superimposed on a much slower oscil- 
lation which has a longer wavelength. The longer oscillation 
is directly linked to the strength of the ferromagnetic compo- 
nent and mainly depends on the singlet part of the anomalous 
Green's function F. The sign of the longer oscillation varies 
with multiple phase transitions from to ir which depend on 
the thickness of the magnetic layer. The rapid oscillations are 
linked to the triplet components f± of the anomalous Green's 
function F. Although of a shorter wavelength to the slower 
oscillations, this rapid oscillation also decays on the scale of 
the magnetic coherence length. 

The main feature of the results presented in this paper 
is that the Josephson coupling through a conical ferromag- 
net may not be long ranged as previously expected. In the 
limit considered, we have shown that the proximity effect is 
short with a length scale comparable to that of the proxim- 
ity effects in a weak and collinear ferromagnet. Thus, the 
theory explained in this paper is complementary to previous 
studie s 36,37 ' 38 which assume that the magnetism of thin-film 
conical magnets is comparable to the magnetism of conical 
magnets in bulk single-crystal form. 

From the experimental view, the theory presented in this 
paper is directly applicable to situations in which two singlet- 
type superconductors are coupled via a rare-earth conical 
magnet (e.g., Ho). The experiment should be designed in 
such a way that the current flowing through the superconduc- 
tor/ferromagnetic/superconductor junction is restricted to the 
growth direction of the conical magnet, e.g., along the z axis, 
such as in the case illustrated in Fig. [TJa). For similar experi- 
mental situations, see references 18 and 49 . 
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APPENDIX: COMPUTATIONAL METHOD IN DETAIL 

The boundary conditions at the S/FM interfaces are given 
by Eqs. (|26| | and d27b . and the equations relating the appropri- 
ate derivatives, 

dz a s dz a s dz dz 

If we assume rigid boundary conditions, dF^ /dz and dF R /dz 
are small because Fl , F R w constant. However, they are 
still not equal to zero and they must not be neglected in 
this treatment. The boundary conditions d261 l. d27l > and dA.lt 
can be written as vector equations: P(0) = M/^F^ and 
F R = M R f¥(d) with the 8x8 matrices 



M 



M Rf = 



Gj 1 ! Rosl 
6 G s a s l/a f 



GA 



Ra f l 



Gj^fl/a, 



(A.2) 



(A.3) 



The symbols 1 and 6 denote the 4x4 unit and zero matrices, 
respectively. In order to deal with the interior of the FM, we 
return to Eqs. ©— dS) and introduce new functions as / + = 
f' + e lQz and /_ = f'_e~ l ® z (primes do not denote derivatives 
here). In this case we obtain the equations 



d 2 fo 
dz 2 



2k 2 f - 2ik 2 J z - ik 2 r {f + + f'_) = 0, (A.4) 



d 2 f! 



^± ± - Q 2 f'± ~ 2*£/£ - 2ik 2 f = 0, (A.5) 



<Pfz_ 

dz 2 



2k 2 J z - 2ik 2 J = 0. 



Equations dA.4t — dA6l > can be written compactly as 

dF'(z) 



dz 



M F F'(z) with M F 



1 
K L 



if we introduce the formal vector 

F'M -ff f f f df ° df + df '- dfz 

The 4x4 matrices K and L are given by 



(A.6) 



(A.7) 



(A.8) 



K 



/2*£ 




iky, Qzkg \ 


2ik 2 r 


Q 2 + 2kl 


1 


2ik 2 r 





Q 2 + 2fcS o 


\ 2ik 2 





2kl ) 



L 





-2iQ 

2iQ 





(A. 10) 



as it is clear from Eqs. dA.4t — dA6i >. Solving Eq. dA.7t in the 
region < z < d gives F'(rf) = exp(dM F )F' (0). We also 
introduce the conversion matrices 





/ 1 o 











\ 








1 


i 

















1 


-i 




















1 













(A. 


f - 








1 












-iC 


1 Q o 





1 


i 








iQ 


Q 





1 


-i 








V 











l) 






/I 




















°\ 





a/2 


1/2q 




















-ia/2 


i/2a 


























1 


























1 














iaQ/2 


-iQ/2a 








a/2 


l/2a 








aQ/2 


Q/2a 








-ia/2 


i/2a 





\o 




















1/ 



Mf F = 



(A.12) 

where a = e l ® d . These can be used to convert F(0) to 
F'(0) as F'(0) = M F/ F(0) and F'(d) to F(d) as F(d) = 
M.fFF'(d). By taking the matrices defined in Eqs. dA.2t . 
(|A31 . dATTT i, (TATTT i and (TA~T2l we obtain 

F R = MF L , (A. 13) 



M = M Rf M fF exp(dM F )M f/ M /L . (A.14) 

Due to rigid boundary conditions, the components of Fl are 
foL = Ae-*/ 2 /VAW and f xL = f yL = f zL = 0, 
while the components of F R are fo R = Ae l *^ 2 /\/A 2 + h 2 u 2 
and f xR = f yR — f zR — 0. The first four components of the 
vectors F^ and F R are therefore known and Eq. dA.13l ) can 
be used to obtain the remaining four components (the deriva- 
tives). If we divide the matrix M into 4x4 blocks as 



M 



Mix M 22 



(A.15) 



and then take the first four components of Eq. dA.13l ), a pre- 
multiplication by gives the derivatives as 



/ df 0L /dz 
df X L/dz 
dfyh/dz 

\ df zL /dz 




( foL 

fxL 

fyL 

\fzL 



(A. 16) 



(A. 9) The first component d/oL / dz is used in Eqs. ( f50b and (fSTb . 
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